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I. INTRODUCTION 


Abstract— A shoaling and breaking water wave model was developed in 
this study using the velocity potential of the solution of the Laplace 
equation. Formulation was carried out using Weighted Kinematic Free 
Surface Boundary Condition Equation, wave number conservation 
equation, and energy conservation equation. In the Weighted Kinematic 
Free Surface Boundary Condition equation, a weighting coefficient was 
used which needs to be determined. The value of the weighting coefficient 
was extracted from the dispersion equation. The dispersion equation was 
formulated using Weighted Kinematic Free Surface Boundary Condition 
and Weighted Momentum Equation to obtain weighting coefficient that 
meets the hydrodynamic equilibrium equation. By using the weighting 
coefficient, good shoaling and breaking results were obtained. 


said that the weighting coefficient is obtained analytically 
and fulfills the law of hydrodynamic equilibrium. 


This study is a follow up research ofHutahaean’s (2021a) 


research. Previously, breaking equations were formulated, Il. 
namely the equations to calculate the breaking wave 
heightH,, breaker wavelength L,and breaker depth hy. 


These breaking equations were 


formulated using 
weighted Kinematic Free Surface Boundary Condition 


POTENTIAL VELOCITY EQUATION 
The total velocity potential equation obtained from 
solving the Laplace equation with Variable Separation 
Method (Dean (1991)) is, 

o(x,z,t) = G(coskx + sin kx) cosh k(h + z) sinot..(1) 


Equation and the weighted momentum equation similar to 


this study. The weighting coefficient in that study was 
obtained by calibrating the breaker heightH, model with 
the breaker height from a number of empirical equations 
of breaker height index. Therefore, the value of the 
weighting coefficient was determined by the breaker 
index equations, instead of based on the hydrodynamic 


equation. 


This study aims to obtain the value of the weighting 
coefficient analytically. The weighting coefficient was 
extracted from the dispersion equation formulated from 
the hydrodynamic equilibrium equations. Thus, it can be 
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x is the horizontal axis;zis the vertical axis; t is time; G 
and k are wave constants, where k is the wave number, 


where L is wavelength k = = o is angular frequency, for 


: 2 
a period of waveT, o = At (1) there are two wave 


constants for which the equation needs to be determined, 
G and k. 

There is a value of kx where coskx = sinkx. This point 
is called the characteristic point on the x-horizontal axis. 
The formulation of the wave constantsGandkis carried out 
at the characteristic point, where the wave constants 
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obtained apply to thecosinecomponent and the sine 
component. At this characteristic point, (1) can be written, 


o(x,z,t) = G cos kx cosh k(h + z) sin ot ...(2) 


It is important to note that there is a double value of G in 
(2). 

I. EQUATIONS OF CONSERVATION 
There are conservation equations in velocity potential (1) 


and (2) (Hutahaean (2020)), namely 


The wave number conservation equation, 


dk(n+3) 

— m 0 s) 

his the water depth, Ais thewave amplitude.In deep water, 

applies 

tank (a +4) = Ch 1 (4) 
2 h esses 


A 
k(h+“) =6n ene (5) 
tanh Om = cp ~ 1....(6) 


0 is a positive number greater than one. The value@is 
determined in the — shoaling-breaking condition, 
where@affects the breaker depth. 


Keeping in mind (3), then (4), (5) and (6) apply to all 
domains, both in deep water and in shallow water. 


The next conservation equation is the energy conservation 
equation, namely, 


ae) i | re (7) 
dx dx 

And 

ag 

0 (8) 


Equation (8) shows that the wave constant G changes 
linearly with respect to the x-horizontal axis. While, the 
water depth changes with respect to the x-horizontal axis, 
thenthe wave constantGalso changes linearly with respect 
to water depth. Furthermore, (7) and (8) indicate that the 
wave numberkalso changes linearly with water depth. 


IV. WAVE AMPLITUDE FUNCTION 


Wave amplitude function is an equation that expresses 
the relationship between G, k and A. This equation is 
formulated using the weighted kinematic free surface 
boundary condition (Hutahaean (2021b)), which takes the 
following form, 
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Y= = Wy - Un ...(9) 


n = (x,t) is water surface elevation equation to the still 
water level; Wy is vertical water surface particle velocity; 


u, is horizontal water surface particle velocity; y is 


n 
weighting coefficient greater than one whose value will 


be determined. 


From (2), the vertical particle velocity is, 


do 


w(x,z,t) = — z —Gk cos kx sinh k(h + z) sin øt 


E (10) 
Vertical water surface particle velocity at z = ņ is, 
W, = —Gk cos kx sinh k(h + n) sin ot AD 
Horizontal particle velocityis 
do 
u(x,z,t) = — ae Gk sin kx cosh k(h + z) sin øt 
sakosi (12) 
Horizontal water surface particle velocity at z = ņ is, 
Uy =Gksinkxcoshk(h+n)sinot se (13) 


Substitution (11) and (13) to (9), 


d 
Va = —Gk cos kx sinh k(h + n) sin ot 


d 
—Gk sin kx cosh k(h +n) sin ot a 


At the characteristic point where cos kx = sin kx 


EE EE EE ce an) 
Ve = —Gk (sinh k(h +) + cosh k(h +1) Z 
cos kx sin ot 

Or, 

LERE E S T EE cake +2) 
Vat = —Gk cosh k(h + n) (tanh k(h + n) += 
coskxsinot > n (14) 


As a periodic function, then, 
d 

Gk coshk(h +n) (tanh k(th+n)+ i = constant 
Thus, (14) can be integrated with respect to time t by 
integrating sin ot only. 

Crs ae Tee V(t hk(h + +5) 
n(x,t) = F cos n) | tan n Ax 
cos kx cos ot 


For 
_ Gk an 
A = “cosh k(h +n) (tanhk(h +n) +2) ..(15) 


The water level equation becomes 
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n(x,t) = A cos kx cos ot ....(16) 


The characteristic point in the time domaintis the value of 
ot where, cosot = sinot. Then, at the characteristic 
point in the x ddomain and in thet domain, 


A 
ae (17) 
an kA 
gee — a (18) 
Substitution (17) and (18) to (15), 
Gk A A kA 
A= —cosh k (h 4 5) (tanh k(n + 5) = <) 
yo 2 2 2 
Substitution (5) and (6), 
ae Gk na ( <) 
= vo T \ Ch z 
Given that there is a double value in G, 
Gk kA 
A = > cosh On (cr -*) tly 


This equation is the wave amplitude function, which can 
be written as the equation for G, 


2yoA 
G = ——— nra 2 
k cosh an(cn—“4) ( 0) 


Ve DISPERSION EQUATION 


Surface momentum equation (Hutahaean (2021b)) 
ignoring convective acceleration, 

2 dün an 
= tet (21) 


Substitutions (13) and (16). The equation is done at the 
characteristic point, 


A 
y7o Gk cosh k (h + 5) = gkA 


kon the left and right sides cancel each other out, and 
thewave amplitude Aon the right side is substituted by 
(19), 


A A Gk kA 
y“oG cosh k (n + 5) = 97, cosh On (cr — < 


The Gon the left and right sides of the equation cancel 
each other out. Keeping in mind (5), the terms of the 
hyperbolic function on the left and right sides also cancel 
each other out, 


2,-— gk _ kA 
y2o= 2 (cn £) sD) 
This equation can be written as, 
Lp- he + y3o2@=0 0 as (23) 
a=% ; b= -7% iS ya" ; d=b?-—4ac 
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—b -vd 
p- 2V1 
2a 


The minus sign used atVd is to make the obtained 
wavelength not too short. 


Table (1) presents the results of the calculation of the 
wavelength with (23), where the wave height H from the 
Wiegel equation (1949,1964)is, 


_ gt? 
~ 15.62 


which is the maximum wave height in a wave period in 
deep water. Assuming a sinusoidal wave, the wave 
amplitude is half of the wave height. 


The calculation was carried out using y = 1.342, the 
and a : H 
determination of y is aimed at getting wave steepness 1 


0.17 which is following the criteria of criticalwave 
steepness from Toffoli et al (2010), where according to 
Toffoli et al, the critical wave steepness can actually 
reach more of 0.2. 


Table (1) Wavelength calculation results using (23) 
T H L 


H 

(sec) (m) (m) L 
6 1.451 8.516 0.17 
7 1.975 11.591 0.17 
8 2.579 15.139 0.17 
9 3.265 19.161 0.17 
10 4.03 23.655 0.17 
11 4.877 28.623 0.17 
12 5.804 34.064 0.17 
13 6.811 39.978 0.17 
14 7.899 46.365 0.17 
15 9.068 53.225 0.17 


Calculation of the wavelength in Table (1) was carried out 
using the value of 0 = 1.95, where the value of 
tanh 6m = 0.999990. However, the value of @ does not 
really affect the wavelength. The effect of 0 is on the 
breaker depth, which is discussed in section (7). 
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VI. DETERMINATION OF THE VALUE 
OF THE WEIGHTING COEFFICIENT 


Y 


Hutahaean (2021b) obtained the value of the weighting 
coefficienty by usingTaylor series on the function in the 
form (1), where the value of the weighting coefficient y = 
3. This coefficient is not yet a product of the 
hydrodynamic equilibrium equation. 


In the previous section, it has been shown that the value 
of the weighting coefficient y can also be obtained by 
adjusting the obtained wavelength with the critical wave 
steepness. However, the value of the weighting 
coefficient is the result of coercion, not naturally 
contained in the hydrodynamic equilibrium law. 


In this section, the value of the weighting coefficient y is 
determined analytically, extracted from the dispersion 
equation (23). This equation is derived from the 
hydrodynamic equilibrium equation. Hence, it can be said 
that it is another form of the hydrodynamic equilibrium 
equation. By extracting the weighting coefficient y from 
(23), the value of the weighting coefficienty is the 
product of the hydrodynamic equation. 


The determinant din(23)must be greater than or equal to 
zero. 


The amplitude of waveAis known variable which 
accordingly taking the equals sign will get the value of 


Ymax: 
2 
Yami =h (25) 


Ifyfrom (25) is used, where the determinant d=0, then the 
wave number equation becomes, 


c 
k=% 
A 


In (26), there is no longer a value of y, but the equation is 
formulated based on the condition of the determinant 
value d = 0, then the value of y in (26) is Ymax. It was 
found that in deep water, the wave-number is only 
determined by the wave amplitude. Meanwhile, the value 
of Cp is always close to or equal to one. Moreover, it is 
also found that the wave number is not determined by the 
wave period T. 


Table (2) presents the results of the calculation of y using 
(25), wave number kis calculated by (26) wave height H 
is calculated by (24), which assuming a sinusoidal wave, 
the wave amplitude A = 0.5H. 
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Table (2) Wavelength is calculated by (26), y is 


calculated by (25) 
T Y H L H 
(sec) (m) (m) L 
6 1.455 1.448 4.55 0.318 
7 1.455 1.971 6.193 0.318 
8 1.455 2.575 8.089 0.318 
9 1.455 3.259 10.237 | 0.318 
10 1.455 4.023 12.639 | 0.318 
11 1.455 4.868 15.293 | 0.318 
12 1.455 5.793 18.2 0.318 
13 1.455 6.799 21.359 | 0.318 
14 1.455 7.885 24.772 | 0.318 
15 1.455 9.052 28.437 | 0.318 


Table (2) shows that value of the weighting coefficient y 
is constant with respect to the wave period of y = 1.455. 
The obtained wavelength is shorter than the previous 
calculation results because a larger weighting coefficient 
is used. the wave steepness obtained is constant for all 


wave periods with a value of = = 0.318, the value of this 


wave steepness is much greater than the critical wave 
steepness value of Toffoli et al. However, the critical 
wave steepness value is the natural value contained in 
(23). 

Equation (26) obtained by working on the weighting 
coefficient ofymaxon (23). Therefore, the calculation of 
the wavelength with (26) at various wave amplitudes 
produce a wave steepness = of 0.318 as presented in 
Table (3). 


Table (3) Wave steepness of (26) 


A H L H 
(m) (m) (m) L 
0.1 0.2 0.628 0.318 
0.2 0.4 1.257 0.318 
0.3 0.6 1.885 0.318 
0.4 0.8 2.513 0.318 
0.5 1 3.142 0.318 
0.6 1.2 3.77 0.318 
0.7 1.4 4.398 0.318 
0.8 1.6 5.027 0.318 
0.9 1.8 5.655 0.318 

1 2 6.283 0.318 
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VII. FORMULATION OF SHOALING- 
BREAKING 
Equation 
The shoaling and breaking equations are formulated using 
the conservation equations contained in (1), (Hutahaean 
(2020)) namely, 


a.The energy conservation equation (7) which can be 


written as 
dG Gdk 

ke ae (27) 
dx 2dx 


b.The wave number conservation equation (3), which can 
be written as 
A\ dk | kdA ah 
(n+ =)S 4525 -K= -...(28) 
2/ dx 2 dx dx 


Equation (19) is differentiated about the x-horizontal axis 
and (27) is substituted for the product of the differential, 


JE (29) 


dx 


dA _ Gcosh6n ( kA 
ax 4yo h 2 


To make writing easier, defined 


_ Gcosh 6n kA 
a = (cn) ....(30) 
Thus, (29) becomes 
dA dk 
an őz tt (31) 


Substituting this equation into (28), it is obtained, 


((n+2)+£a)%=-42 32) 


dx dx 
: PARR? ak ak. ; 
This equation is used to calculate ax After ax IS obtained, 
aA : 
gem be calculated using (29). 


Calculation steps: 


1.At a pointx = x, givenh,, Ay, Ky, Gx, 
will calculate Ay45,, Kys5,and Gy+5x, at pointx = x + ôx 


-o — GxB _ kxAx 
2.Calculate :@ = rice (cn > ) 


3. Calculate = with the following equation : 


Ee dk __„ ah 
(rx 2/724 jax “ax 
dk 


aA. . aA 
4. Calculate —with the equation: — = a — 
dx dx dx 


5. Calculate k,.5,andAy45, : 
dk 
a. Ky 5x = kx + Ôx 
aA 
b. Ayssx = A, + 6x 
6. Gy+sx1S calculated analytically by integrating (27), 


1 
G tox = en Gy—F (Ink x4 8x —Inkx) 
x+6x 
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At the starting point, namely in deep water, wave 
numberkwas calculated by (26), Gwas calculated by (20) 
and wave amplitudeAwas calculated by (24). 


Fig (1) visualizes the results of the shoaling-breaking 
model for a wave period of 8 sec., the deep water wave 
height Hp was calculated by (24) assuming a sinusoidal 
wave was A= 0.5Họ. As a calculation parameter 
ofy(25) and the coefficient ofdeep water depth 0 = 1.95 
were used where tanh 07 = 0.999904.The results of the 
calculation of H, are presented in Table (3), namely Hp = 
2.987 m while the results of the calculation using (33), 
namely the equation of Komar&Gaughan (1972), 


obtained H, = 3.019 m. Breaker depth index ae = 0.78, 
b 


according to Mc Cowan’s (1984) breaker depth index 
equation. 


4 
3 
£2 
je By 
1 
0 
0 2 4 6 8 
h (m) 


Fig (1) Example of a shoaling-breaking model result 


The equation for breaking wave height H, from Komar 
& Gaughan (1972) is as follows: 


H, = 0.39 g 5 (TRY 5 eee (33) 


This equation uses a primitive variable, namely the wave 
period T, making it very practical in its use. 


0.04 
T 0.03 
fd) 
< 
N 0.02 
E£ 
© 0.01 
0 
0 2 4 6 8 
h (m) 


Fig.(2). Change ofGwith water depth h 


(8) shows that the change in the wave constant Gwith 
respect to water depth is linear. In Fig (2), the linear 
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nature only occurs before breaking, while after breaking, 
it is nonlinear, especially in very shallow waters. 


10 
8 
= 6 
E 
2 4 
2 
0 
0 2 4 6 8 
h (m) 


Fig.(3). Change of wavelength Lwith water depth h 


Using (7) and (8), it was found that the change in wave 
numberkon water depth is linear, the same with the one 
on wavelengthL. Figure (3) shows that the change in 
wavelengthLis linear, but slightly nonlinear in very 
shallow waters.With changes in G and changes in k as 
shown, further research and development is needed on 
wave conditions after breaking, especially in very shallow 
waters. 
Table (4) Results of the shoaling-breaking model for a 
number of wave periods 


T Ho H, (m) Hp Hp 
(sec) | (m) | model K-G Ly h, 
6 1.451 | 1.681 1.698 | 0.637 | 0.781 
7 1.975 | 2.287 | 2.312 | 0.638 | 0.782 
8 2.58 | 2.987 | 3.019 | 0.637 | 0.781 
9 3.265 | 3.78 3.821 | 0.637 0.78 
10 | 4.031 | 4.667 | 4.718 | 0.637 | 0.781 
11 | 4.878 | 5.646 | 5.708 | 0.637 | 0.781 
12 | 5.805 | 6.719 | 6.794 | 0.637 | 0.781 
13 | 6.812 | 7.886 | 7.973 | 0.637 | 0.781 
14 | 7.901 | 9.146 | 9.247 | 0.637 | 0.781 
15 9.07 | 10.499 | 10.615 | 0.637 | 0.781 
Note : K-G, Komar&Gaughan 


The stability of breaking characteristics including 2 and 
b 


2 in the model on the wave period is shown in Table (4), 
b 


ie Rie fo H H : 
indicating that both and oa can be said to be 
b b 


constant.The calculations in Table (4) were carried out 
usingyof (25), the coefficient ofdeep water depth 0 = 
1.95. 


www.ijaers.com 


International Journal of Advanced Engineering Research and Science, 9(6)-2022 


To determine the effect of the deep water depth 
coefficient@, a calculation was carried out using 0 = 1.85 
with y from (25). 


Table (5) Results of the shoaling-breaking model with 
0 = 1.85 


T Ho H, (m) H, H, 
(sec) | (m) | model | K-G Ly h, 


6 1.451 | 1.681 | 1.698 | 0.639 | 0.835 


1.975 | 2.287 | 2.312 | 0.638 | 0.833 


2.58 | 2.987 | 3.019 | 0.637 | 0.832 


S| coo}; N 


3.265 | 3.78 3.821 | 0.637 | 0.833 


10 | 4.031 | 4.667 | 4.718 | 0.638 | 0.833 


11 | 4.878 | 5.646 | 5.708 | 0.637 | 0.831 


12 | 5.805 | 6.72 6.794 | 0.637 | 0.832 


13 | 6.812 | 7.886 | 7.973 | 0.637 | 0.832 


14 | 7.901 | 9.146 | 9.247 | 0.637 | 0.831 


15 9.07 | 10.499 | 10.615 | 0.637 | 0.832 


The calculation results in Table (5) present that there is no 
change in the breaking wave heightH,, even though the 
deep water depth coefficient @was reduced. The changes 


H H : 
that occur were the enlargement of z and a This shows 
b b 


that the reduction of the deep water depth 
coefficientOcauses the shallower breaker depth hy. 


VIII. CONCLUSIONS 
This study concludes that the weighting coefficientyon 
theweighted kinematic free surface boundary condition 
and on the weighted momentum equation can be extracted 
from the dispersion equation. 
The dispersion equation is formulated using 
hydrodynamic equilibrium equations. Accordingly, the 
value of the weighting coefficient obtained is to meet the 
laws of hydrodynamic equilibrium equations. 
A good shoaling-breaking model is obtained using the 
weighting coefficient extracted from the dispersion 
equation. However, a large critical wave steepness is 
obtained both in deep waters and at the breaking point. 
Therefore, further research is still needed, both on the 
model and on the existing critical wave steepness criteria, 
both in deep water and at the breaking point. 
In general, the results of the shoaling-breaking model 
provide good information regarding the height of the 
shoaling wave and the breaking wave height. However, 
further research is still needed on wave conditions after 
breaking, especially in very shallow waters. 
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